Collectively-enhanced optomechanical coupling in periodic arrays of scatterers 
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In a recent publication [A. Xuereb, C. Genes, and A. Dantan, Phys. Rev. Lett. 109, 223601 
(2012)] the optomechanical properties of a periodic array of identical scatterers placed inside an 
optical cavity were investigated. It was shown that operating at a point where the array is trans- 
missive results in a linear optomechanical coupling strength between the cavity field and a collective 
motional mode that may be several orders of magnitude stronger than is possible with an equivalent 
reflective ensemble. In this publication we describe and interpret these effects in significantly greater 
detail, investigate the nature of the scaling laws noted previously, and discuss their robustness to 
imperfections. 
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The ability to measure and control the motion of mas- 
sive mechanical oscillators has progressed dramatically 
in recent years [1, 2], and several important milestones 
have been reached in the field of optomechanics towards 
bringing this level into the quantum regime, including the 
cooling to the motional quantum ground state [3-5] , the 
detection of quantized mechanical motion [6, 7], and the 
observation of the ponderomotive squeezing of light [8, 9] 
or of the radiation-pressure shot noise on a mechanical 
oscillator [10]. 

One challenge faced by the current generation of op- 
tomechanical experiments is that the interaction strength 
between a single photon and a single massive mechani- 
cal element is typically very weak. This can be ame- 
liorated by confining light in wavelength-scale struc- 
tures [11] or, generically, counteracted by the use of 
strong light fields in an optical resonator to amplify the 
interaction strength [12], albeit at the expense of trad- 
ing off the intrinsically nonlinear nature of the radiation- 
pressure interaction (see, however, the recent proposals 
in Refs. [13, 14]). A growing number of theoretical pro- 
posals has contemplated the opposite - 'strong coupling' 
- regime, where a single photon can affect the motion of 
the oscillator significantly, thus giving access to the full 
quantum nature of the optomechanical interaction [15- 
23]. 

On the other hand, collective effects in optomechan- 
ical systems involving multiple mechanical and electro- 
magnetic field modes have been discussed in a number of 
theoretical works, in connection with, e.g., optomechan- 
ical entanglement [24-34], enhanced displacement sensi- 
tivity [35-37], optomechanical nonlinearities [20, 21, 38- 
40], quantum information processing [41-44], many-body 
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physics [45-49] , as well as in a number of recent experi- 
ments [50-56] . Collective optomechanical effects are also 
at the heart of cavity optomechanics with cold atomic 
ensembles [57]. 

Motivated by the exploration of such collective op- 
tomechanical effects, we recently [58] showed that the 
collective motion of a periodic array of identical scatter- 
ers, when placed inside a cavity field, can couple very 
strongly to the optical field in the configuration where 
the array is transmissive, in contrast to the usual reflec- 
tive optomechanics approach. The aim of the present 
work is to present a more detailed exploration of this 
system, to highlight the regimes in which these generic 
collective effects are seen, and to explore its robustness 
to imperfections. 

This paper is organized as follows. In the next section 
we discuss the tools used to model a periodic array of N 
identical elements, and show that such an array can be 
modeled as a single effective element in the framework of 
the transfer matrix theory for one-dimensional scatter- 
ers. Section II discusses the optomechanical properties 
of such a generic iV-element stack, when placed inside 
an optical cavity, in two distinct and opposite regimes: 
(i) A maximally-reflective stack (Sec. II B), and (ii) A 
transmissive stack (Sec. TIC). We also show that the 
equations of motion for the system at hand are effectively 
described by an optical field interacting with a single col- 
lective mechanical mode. The next section discusses the 
second regime in further detail, and explores the scal- 
ing of the optomechanical coupling strength in the trans- 
missive regime with the various system parameters; we 
show that the increase in the optomechanical coupling 
strength with the number of scatterers combines with an 
effect whereby the linewidth of the cavity resonance is 
narrowed due to the presence of the stack to provide an 
enhancement of the optomechanical cooperativity by sev- 
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FIG. 1. Schematic diagram illustrating the basis of the model: 
A^ scatterers interacting with four running waves. B and C 
represent incident fields, A and D outgoing fields. The four 
field amplitudes are connected by means of transfer matrices. 



eral orders of magnitude over that of a single element sys- 
tem. Section IV examines the robustness of these effects 
with respect to deviations from the ideal periodic array 
of identical, lossless scatterers, and Sec. V provides nu- 
merical estimates for the mechanisms we predict, based 
on state-of-the-art experiments. 



I. MODELING AN iV-ELEMENT STACK 

Throughout this paper we shall restrict ourselves to 
a one-dimensional system and use the transfer matrix 
formalism to model a periodic TV-element array, as illus- 
trated in Fig. 1. As is well-known [59, 60], this formal- 
ism can treat elements ('scatterers') that interact linearly 
with the electromagnetic field, no matter the strength of 
this interaction or whether the scatterer is lossless or not. 
Within this formalism, each scatterer is parametrized by 
its polarizability C, which is real for lossless scatterers 
but complex in the presence of absorption. Q is related 
to the amplitude reflectivity r of the element through the 
relation 



< 



l-< 



and, for a lossless scatterer, we have 



i+e 



(1) 



(2) 



which allows us to link our results to, e.g., those in 
Ref. [61] through their Eq. (8). The flexibility afforded by 
the transfer matrix formalism allows us to treat ensem- 
bles of atoms in an optical lattice on the same footing as 
periodic arrays of macroscopic scattering elements (e.g., 
arrays of thin dielectric membranes) in the limit of a one- 
dimensional scattering theory, as we illustrate in App. A 
and App. B. 



A. Transfer matrix for A'^-element stack 

Let us consider N equally-spaced, identical and non- 
absorbing elements, each of which has a thickness much 
smaller than an optical wavelength. The spacing d be- 
tween these elements determines the overall optical prop- 
erties of the ensemble. Our starting point is the matrix 



that links the fields interacting with a single element of 
polarizability C, [59], 



A/m(C) = 






(3) 



and the matrix that describes propagation of a 
monochromatic beam of wavenumber k over a distance d 
through free space. 



Mp(d) 



„ikd 









— ikd 



(4) 



both of which have unit determinant. These matrices re- 
late forward- and backward-propagating electromagnetic 
waves on either side of the element: 



M 



(5) 



where A and C {B and D) are the amplitudes of the 
backward- (forward-) propagating waves, as illustrated 
in Fig. 1. We wish to evaluate a product of the form 



M^(C)-Afp(d)-Af,n(C)---M„,(C), 
with N copies of Af,n(C)- First, we note that 

Mp{d/2) ■ M^iC) ■ Mp{d) ■ ■ ■ M^iC) ■ Mp(d/2) 

= [Mp(d/2) • Af,„(C) • Mp{d/2)f = M^ 

where the second line defines the matrix M: 



(6) 



M 



(l + iQe''"^ iC, 

-iC (1 - iC)e-''"' 



(7) 



(8) 



We can easily see that det M — 1, whereby it can be 
shown [62] that for real C we can write 



M^ = 



(1 + ix)e'('"^+t''> 
-iX 



(1 - ix)e-''-'"^+'''^ 



(9) 



where x = Ct^w-i(a), with C/„(x) being the n^^ Cheby- 
shev polynomial of the second kind, a — cos{kd) — 
Csin(/cd), and 



1 - iCUN-iia) 



(l-zC)t/jv-i(a)-e*'=^C/jv-2(a) 



(10) 



Upon removing the padding of d/2 from either side, we 
obtain the matrix that describes the iV-element ensem- 
ble: 



Mn = Mp[fi/{2k)] ■ M,n(x) • Afp[M/(2fc)] 



(11) 



In this respect, N lossless elements behave as a collective 
'superelement' of polarizability x along with a phase shift 
/i/2 on either side of the stack. 
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FIG. 2. Intensity reflectivity of A'^ = 5 lossless elements, each 
of which has reflectivity 20% (shown as the dotted green line), 
as the spacing between the elements is scanned. We denote 
the inter-element distance where the reflectivity attains its 
peak as do, and two of the distances where the reflectivity is 
zero as d±. The curve is periodic with period A/2. 



Let us note, finally, the link between the transfer ma- 
trix of any optical system and its amplitude transmissiv- 
ity and reflectivity. Indeed, suppose we can describe a 
system by means of the transfer matrix 



mil ™i2 

m21 TO22 



(12) 



Then, the complex transmissivity of the system is simply 



T 



1 

TO22 



and the corresponding reflectivity 



n 



■mi2 
■m22 



(13) 



(14) 



These quantities are used throughout this paper to char- 
acterize the optical properties of our system. A plot of 
the intensity reflectivity of a lossless 5-element ensem- 
ble as the spacing between pairs of elements is varied is 
shown in Fig. 2. Note that, despite each element having 
a reflectivity of 20%, the reflectivity of the entire ensem- 
ble ranges from 0% (e.g., at d = d±) to about 97% (at 
d^do). 



II. OPTOMECHANICS OF A PERIODIC 
ARRAY OF SCATTERERS 

Our next task is to place the array just described in- 
side a near-resonant optical cavity. The interaction of the 
array with the cavity held will shift the resonances of the 
main cavity. As is usual in optomechanics, we assume 
that each element in the array is harmonically trapped, 
and are interested in one particular figure of merit: the 
coupling strength equal to the frequency shift incurred 
by the cavity resonance when the array undergoes a dis- 
placement equal to the size of the relevant zero-point- 
fluctuations. As a yardstick we shall use the quantity 



2u}rX 



9 



zpt 



L 



(15) 



which is the optomechanical coupling strength for a per- 
fectly reflective mirror near the center of a cavity of 
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FIG. 3. Schematic of the model we consider. A'^ equidis- 
tant elements, each independently harmonically-bound, are 
positioned in the field of a Fabry-Perot resonator (we shall 
consider only the case for which L ^ Nd). 



length L and resonant frequency lOc [63]. The size of 
the zero-point-fluctuations of each element is denoted by 
Xzpt = \/fi,/(mWm), where m is the effective mass of the 
element and Wm its oscillation frequency. 



A. Element stack inside a cavity 

Figure 3 illustrates schematically the periodic element 
array placed inside a Fabry-Perot cavity of length L, 
assumed much longer than the array (L ^ Nd). The 
transfer matrix describing this system is then 



Meav = M„,(Z) • Mp(L/2 + x) 

■Mn • A/p(L/2 



x)-M^iZ). (16) 



Here x is the displacement of the ensemble with respect 
to its position at the center of the cavity and Z is the po- 
larizability of the cavity mirrors, assumed equal for both. 
For good, lossless, cavity mirrors {\Z\ > 1), the finesse of 
the cavity may be simply written J- — Tr\Z\\/Z^+l [64]. 
The transmission of the system, following Eq. (13), is 
given by 



/cav — 



1 



Mr: 



(17) 



the maxima of T^av give the resonances of this system. In 
order to find these resonances analytically, we consider a 
simpler system where the cavity mirrors are perfect; we 
need only solve the relation 



1 
-1 
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e-'^_ 



1 
-1 



with 9 = k{L/2 + x 
We thus obtain 



(18) 
+ x) + /i/2 and = k{L/2 - x) + /i/2. 



AkL 



-tfj. 



i + «x 



iXC0s(2fca;) ± V 1 + x^ sin^ C^kx) 



(19) 



However, we immediately see that this equation is tran- 
scendental in fc, and therefore cannot be solved analyti- 
cally; this equation is easily solvable for L, however, and 
it therefore easy to solve this equation to find the reso- 
nant length of our cavity. 
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FIG. 4. Transmission through a cavity with N — 5 immo- 
bile elements configured for maximal reflectivity. The dashed 
lines denote the bare-cavity resonances, which are shifted due 
to the presence of the ensemble, x is normalized by a factor 
yN, such that the gradient of the bright curves gives directly 
the linear optomechanical coupling at that point. ((" — —0.5, 
L ~ 6.3 X lO'^A, d = d-, bare-cavity finesse « 3 x 10'*, corre- 
sponding to cavity-mirror reflectivities of 99.99%.) 



where 7^ = x§/(l + Xo)- The factor of 1/VN that 
appears in this expression has a natural explanation: 
The motional mass of N elements is Nm, and there- 
fore the scale of the zero-point- fluctuations is x^pt/VN- 
We can draw two immediate conclusions regarding gcom'- 
(i) 5com < g, and (ii) ^com is optimized for 7?. — > 1. A plot 
of the transmission spectrum of a cavity with a 5-element 
ensemble inside it is shown in Fig. 4 as the position of the 
ensemble is varied. The gradient of the bright curves at 
each point gives directly the linear optomechanical cou- 
pling strength at that point. 

The essential aim of this paper is to outline a mecha- 
nism [58] whereby coupling strengths much larger than g 
can be obtained, despite keeping Wc and L fixed. To do 
this, we shall now explore the coupling strength to the 
motion of individual elements, rather than to the ensem- 
ble as a whole. 



B. Center-of-mass coupling: Reflective 
optomechanics 

It is now a legitimate question to ask: 'If d (or x) shifts 
by a small amount, how much will the resonant frequency 
of this cavity shift?' This question is, of course, easily an- 
swered by expanding Eq. (19) in small increments about 
its solution. Assuming a dominantly linear effect, we 
replace fc ^ fco + Sk, x ^ x + Sx, x ~^ X + ^X^ and 
fi ^ fi + 6fi in Eq. (19). Around resonance, the result 
simplifies to 



L6k + Sfi = 



-1 ± cos{2kox)/^Jl + x^ sin^(2fco2:) 
x5x/(l + x') 



T 



2x sin(2fcox) / y 1 + x^ sin^ (2fcoa;) 
X [x5k -|- ko6x) . 



(20) 



We can now use Eq. (20) to compute the coupling 
strength for the center-of-mass motion of the ensemble. 
For such a uniform displacement, d^j, = dx = 0, and we 
assume that \L/x\ is very large, such that we can write 



LSk^T 



2xsin{2kQx) I \J I + x^ &\T?{2kQx) 



ko6x . 



(21) 

The right-hand-side of this equation is maximized when 
sin(2fco2;) — Tlj whereby 



LSk = 2ko{-x/\/lTx^)Sx. 



(22) 



This is, in absolute value, a monotonically-increasing 
function of |x| and is therefore maximized when x 
attains its largest value, xo = CUn-i{\/^ + C^) = 
—i sin[A^ cos^^ (■\/T+^)] . The ensemble attains this re- 
flectivity for when kd = kdo = — tan~-'^(C) (see Fig. 2). 
The resulting coupling strength can then be shown to be 

gcom = gVn/N, (23) 



C. Coupling to each individual element in the 
transmissive regime 

Just as the ensemble attains its peak reflectivity for 
d = do, we can see from Fig. 2 that its reflectiv- 
ity is zero at d = d±, where kd± = — tan^^(C) ± 
cos"-^ [(1 -I- C^)^"^^^cos(7r/A^)]. We now work with this 
inter-element separation and obtain the optomechanical 
coupling strength, i.e., the shift in cavity resonance fre- 
quency, due to the motion of each element in the en- 
semble. To allow one element, say the j* , to move in- 
dependently of the rest of the ensemble, we conceptually 
split the ensemble into three sections: the elements to the 
'left' of j, the j**^ element itself, and the elements to the 
'right' of J. With this logic, the matrix Mjy representing 
the ensemble can be written, for 1 < j < N , 



1 + < < 



*^i/2 1 \l + ixi iXi 
e"*^i/2j [ -ixi 1 - iXi 

^-i{fj,i/2+v+kSxj) 

'^i{fj.2/2+iy-kSxj) Q 

e 

gV2/2 

e~*^=^/2 



1 + iX2 iX2 
-iX2 1 - iX2 



i{ti2/2+iy — kSxj) 

(24) 



where v = kd±, ni and xi describe the ensemble formed 
by the rii = j — 1 membranes to the 'left' of the j*^, and 
/i2 and X2 the one formed by the 112 — N — j membranes 
to its 'right'. The displacement of the j"^ element is 
denoted 5xj ; all other membranes are in their equilibrium 
position. This small displacement shifts the resonance 
frequency of the cavity w — >■ w — gjSxj, defining gj as the 
optomechanical coupling strength for the j"^ element. In 
the transmissive regime, to lowest order in kSxj in each 
entry, the matrix product above can be written, with the 



M 


0.5 - 




- 


<3 


-0.5 - 




-1 - 



:■■""*■■. 

*.. X :. 



2 3 4 5 

Element number, j 



I 




FIG. 5. (Top) Individual coupling strengths for the case of 
5 elements; the dashed curve is drawn as a guide to the eye. 
(Bottom) An illustration of the motional mode that is excited 
in the ensemble, corresponding to the coupling strengths in 
the upper figure. We use membranes to make our illustration 
concrete here, but all our considerations hold for any periodic 
array of scatterers. 



above choice for u. 



1 5xi 



(3 5xj 



(3* Sxj 



-i^ 



a 



6xi 



(25) 



where a and /3 are increments of first order in the relevant 
displacement [note that the (ofF-)diagonal terms are com- 
plex conjugates of each other; this is different to the case 
where absorption is nonzero] . When this matrix is substi- 
tuted into the equation for the resonance condition, the 
terms involving Re{e~'^a| and Re{/3} drop out entirely 
for a symmetric system, such that it suffices to consider 
only the imaginary part of the increment. Let us reiter- 
ate that this happens only because the off-diagonal terms 
are complex conjugates of each other; were absorption to 
be nonzero, this would no longer be the case. Eq. (20) 
now simplifies to 



dk 

86 Xj 



fm{^- 



4 



2d|i 

ov 



(26) 



with 



a = 2ikC 



e*''Kl + iXi)X2-e*'^^Xi(l + %) , (27) 



and 



/3 = 2fccfxiX2 - (1 + *Xi)(l - iX2)e^^^^"^^^l , (28) 
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Sinusoidal mode displacement, x [units of A] 

FIG. 6. Transmission through a cavity with N = 5 immo- 
bile elements in the transmissive regime. This figure should 
be compared to Fig. 4. x is normalized by a factor yjN/2. 
(Parameters as in Fig. 4.) 



which we can rewrite, by expressing xi,2 and /ii.2 in terms 
of Chebyshev polynomials, as 



a = 2ikC 



{I - iQUn,-i{a) - e'''Un,-2{a) 

(1 - iC)C/„3-i(a) - ei^Un^-2{.a) 



(29) 



and 



/3 = 2kC\ C'[/„,_i(a)[/„,_i(a) - [l + C'C/',-i(a)] 

^ {l-iC,)Un,-i{a)-e'''Un,-2{a) 
"" {1 - iOUn,-i(a) - e'''Un,-2{a) 



(30) 



These two expressions simplify considerably to yield 



Iml/S + e-'^a} =2fcCcsc(^ 



iv) 



sm-|-)+C2-C 



sin 27r 



N 



(31) 



and therefore, when C < 0, 



9j 



Ccsc(^) 



-2uicXq 

L - 2Nd( csc^ 






X sin 



, N 



(32) 



We illustrate the 'profile' of the gj foiN^d in Fig. 5. 
As we shall discuss below, the coupling of the collective 
motion of the membranes to the cavity field is governed 

by the constant gsin = \/ Si=i ff?' such that for N — 2 



ffsin 



V2g<:{^/T+e-c) 



(33) 



and for A^ > 2 

fW 5Ccsc(^) 



'-i7r)+e'C 



l-2iVfCcsc2(f)Jsin2(f)+C 



(34) 



This two expressions do not agree upon setting A^ = 2 in 
the latter; this anomaly is due to the relation 



\ 



N 



E^ 



sm^ f 27r'^ 



N 



for TV = 2 
for iV > 2 



(35) 



Plots of the type of Fig. 6, similarly to Fig. 4 but now 
for the case of the sinusoidal mode, can be used to nu- 
merically extract gsim by measuring the gradient of the 
bright curves. Such data agree with the analytic results 
derived above. 



D. Collective motional- mode treatment: 
Heisenberg— Langevin formalism 

In the preceding section we derived a coupling con- 
stant gj that relates the motion of the j"^ element to the 
resonance frequency of the cavity. A particular feature 
of multi-element arrays is that the cavity field couples 
to a collective motion of the elements, with the gj play- 
ing the role of choosing the 'profile' of the mode that is 
coupled to the cavity field, in the spirit of Fig. 5. To 
see this, let us describe the motion of the j"^ mechanical 
element (1 < j < N) through the annihilation opera- 
tor bj , which obeys the Heisenberg-Langevin equation of 
motion [65, 66] 



bj = -{iujrn + T)bj + Fj + V2rfj- 



(36) 



where ^j is the relevant Langevin noise term whose prop- 
erties we leave unspecified. For simplicity, we assume 
that all the oscillators have identical oscillation frequency 
uim, decay rate F, and temperature T, such that in ther- 
mal equilibrium they all have the same average occupa- 
tion. Fj = gjF is a force term due to the action of the 
cavity, whose exact form is not relevant here. To describe 
the collective motion, we use the vector (gj = gj/gsin), 
which is naturally normalized such that X]i=iff? ~ 1' 
and define: b = J2i=i 9j^j ^^'^ ^ — J2j=i 9j^j- Thus: 



P- 



-(iw„i + r)6 + g,i„F + V2re 



(37) 



Under the assumption that the noise terms ^j are of a 
similar nature to one another and are independent (i.e., 
any cross-correlator between ^j and ^j is zero for i ^ j), 
then ^ obeys the same correlation functions as each indi- 
vidual noise term, because of the normalization of (gj), 
whereupon b behaves as a single collective oscillator with 
decay rate F. Let us remark at this point that our de- 
scription in terms of this collective mode is one where we 



merely rotate to a different basis for this A^-dimensional 
space, and therefore the correct normalization, necessary 
for the rotation to be a unitary operation, is J2j=i d] ~ ^■ 
Therefore, the dynamics of the cavity-mechanical sys- 
tem can be described entirely through an optomechanical 
Hamiltonian connecting a single cavity mode to a single 
collective mechanical mode with coupling strength gsim 
mechanical frequency uj^, decay rate F, and noise oper- 
ator ^. 

In the next section we investigate some properties of this 
coupling strength gsin, and draw general conclusions re- 
garding the optomechanical coupling of a transmissive 
ensemble. 



III. TRANSMISSIVE OPTOMECHANICS 

Thus far, we have derived an expression for the cou- 
pling strength of a cavity field to a periodic array of ele- 
ments whose compound reflectivity is zero at the cavity 
frequency. In this section we shall discuss the result- 
ing collective optomechanical coupling. We shall then 
proceed to discuss one further important consequence 
of this collective coupling: a linewidth-narrowing effect 
where the effective linewidth of the cavity resonance de- 
creases over its 'bare' value as the number of elements, 
or their reflectivity, is increased. 



A. Enhanced optomechanical coupling 

The first thing we note is the fact that ^sin is no longer 
bounded above by g. Indeed, for large A^ and |C| but 
small d/L, the expression (34) for gain simplifies consid- 
erably to yield 



n/2 



5^2^3/2 



l + ^K'N^ 



n 



(38) 



We shall provide numerical examples later to show how 
dsin/g can be of the order of 10^-10'*. To explore the scal- 
ing of gsin/ g with A^ we begin by considering a long cavity 
(L 3> d) and approximate the denominator of Eq. (34) 
by 1, obtaining 



gsin = ^/N/29\<:\ [V^TWUW - mM] : (39) 

for large A^. One can now distinguish between two cases: 
(i) {N\C\/n) « 1: gsin/g = \C\VW^^ yielding the VN- 
scaling observed in, e.g., atom-cavity optomechanics ex- 
periments [67-70] that involve large ensembles of very 
low reflectivity scatterers. The coupling strength is also 
multiplied by a factor of |C|, which amplifies the inter- 
action when 1^1 > 1; both these features represent a 
markedly different behavior from the reflective regime. 
(ii) (A^ICI/tt) > 1, which gives the A^3/2_ ^^^^^ (2_g(,g^jjj^g 
shown in Eq. (38). This scaling with N^''^ is a conse- 
quence of the modification of the field mode profile in- 
side the cavity. Indeed, as A^ increases the fraction of 
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FIG. 7. Collective coupling strength gain (blue circles, closed 
red squares, orange triangles) for a cavity of length L « 
6.3 X lO'^A and several choices of scatterer reflectivity, com- 
pared to the center-of-mass coupling strength gcom (open red 
squares) and g (green line). For the red, orange, and blue 
data points we choose a per-element intensity reflectivity of 
20%, 50%, and 99.4%, respectively. For the center-of-mass 
data we illustrate the N~^'^ scaling law that applies for large 
A'^, whereas for the sinusoidal mode we draw a curve through 
the data points as a guide to the eye. Throughout this plot 
we take d — d-, except for the data denoted by the open blue 
circles, for which d = d- + 20A. 
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FIG. 8. (a) Optimal number of elements A'^opt, as a function 
of the single-element reflectivity. For reflectivities close to 
100%, this number decreases rather quickly to < 10. (b) The 
sinusoidal coupling, normalized to g, optimized as A^ = A'^opt 
(solid blue curve), compared to the case for N = 2 (dashed- 
dotted red) and Pcom (dotted orange). For these plots we used 
L^6.3x 10*A andd = d_. 



the energy density per photon inside the array increases 
strongly; an increased optomechanical coupling strength 
is consistent with this increase of energy density [71]. The 
former case, on the other hand, corresponds to an essen- 
tially unperturbed cavity field mode, where there is no 
such concentration of energy density, and a weaker scal- 
ing with N is therefore observed. 

Similarly, in the behavior of gsin/g as a function ^ for 
N = 2 and in the same 'long cavity' limit, one can dis- 
tinguish between two regimes [cf. Eq. (33), setting the 
denominator to 1]. For |^| <^ 1, gsin/g grows linearly 
with Id as a consequence of the increased reflectivity of 
each element. For \(\ ^ 1, however, the cavity field mode 
is modified substantially and is strongly concentrated in 
the region between the two elements. This concentra- 
tion grows quickly as a function of |C| and gives rise to a 
quadratic scaling oi gain/g with |C|. 

The denominator in the full form of Eq. (38) can be inter- 
preted as an effective renormalization of the cavity length 
from L to Lcs ^ L + ^dC/^N^. We shall discuss the 
regime in which the effective length is significantly larger 
than L, where an interesting linewidth-narrowing effect 
occurs, in greater detail in Sec. IIIB below. We have 
already seen that for N , |C|, and d/L small enough that 
Lgff « L the coupling strength scales as (^'^N^^^. On the 
other hand, when the parameters are such that L^s ^ L, 
gsin does not depend on ^ and decreases as N~^/^. We 
say that, in transitioning between the two scaling laws, 
5sin saturates (i.e., reaches a maximum value at some fi- 
nite value for N) before it starts decreasing. 



Optimizing gsin over TV for arbitrary L/d, in this man- 
ner, we obtain g°P^ = ^g^yL/d\(^\. This expression is 
valid for \(\ that is not too large, since the optimal num- 
ber of elements must be > 2. This favorable scaling with 
both N and |C| is a significant improvement over the state 
of the art. Close inspection reveals that g°P^ is propor- 
tional to l/yLd and therefore can be improved either by 
making the main cavity smaller (i.e., decreasing L) or, in- 
dependently, by positioning the elements closer together 
(decreasing d to a smaller value whilst maintaining the 
condition of zero reflectivity). 

Alternatively, one may optimize the parameters such that 
the second term in the denominator of Eq. (38) domi- 
nates, i.e., where icff S> L. The coupling strength then 
takes the approximate form 



V2-'d 



g^N-^^' 



N=2 



(40) 



where we have taken the optimal {N = 2) case and de- 
fined 5sm = i^cXzpt/d as the optomechanical coupling 
strength for a (small) cavity of length d with a single 
moving mirror. In this regime, therefore, the system acts 
as a small cavity of length d and is sensitive to relative 
motion between the two elements but not to the length 
of the main cavity. Figure 7 illustrates, primarily, the 
enhancement of optomechanical coupling strength that 
can be obtained by operating in the transmissive regime 
as compared to coupling to the center-of-mass motion. 
Several important observations can be made from this 
figure. The largest value for ^sin increases with |C|; in 



the figure we show data points for C — —0.5 (20% inten- 
sity reflectivity; red data), C = —0.5 (50%; orange) and 
C — —12.9 (99.4%; blue), the first two of which represents 
a typical reflectivity for SiN membranes used for optome- 
chanical experiments [72], and the last membranes with 
increased reflectivity due to the use of sub-wavelength 
patterning [73, 74]. Secondly, the value of N for which 
the coupling strength is optimized is highly dependent 
on the value of d\ the reflectivity of the array depends 
on d mod (A/2), so that one is free to increase the el- 
ement spacing by integer multiples of half a wavelength 
without affecting its transmission properties (in the limit 
of a wavelength-independent reflectivity). However, the 
coupling strength is sensitive to this increase: The larger 
djL is, the earlier the saturation point is reached, beyond 
which increasing N lowers the coupling strength. The 
figure also illustrates two scaling laws that we derived 
above. The center-of-mass coupling decreases, for iV > 3, 
as N^-^^^, whereas the sinusoidal coupling strength ^sin 
scales approximately as N^/^ for N large enough, but also 
small enough to avoid the effects of saturation. These 
two scaling laws are illustrated by the red dotted curves, 
drawn as guides to the eye. By contrast, the other curves 
and all the data points are generated using the full an- 
alytical formulae, which are in excellent agreement with 
numeric calculations. 

In Fig. 8 we study (a) the optimal number of elements 
required, iVopt, and (b) the resulting coupling, as a func- 
tion of the polarizability of each element in the array. 
For weakly reflective elements, as illustrated in Fig. 8(a), 
the coupling only saturates for very large values of N, 
whereas as the reflectivity of the elements increases the 
optimal number of elements decreases, at first steadily 
and then quite sharply, until it reaches a point where the 
coupling strength decreases for A^ > 2. This curve is sen- 
sitive to the ratio d/L; smaller values of this ratio result 
in larger values for iVopt,. In panel (b), we illustrate the 
optimized sinusoidal coupling strength (solid blue curve), 
as well as the coupling strength for A^ = 2 (dashed-dotted 
red) and the center-of-mass coupling (dotted orange). 
Two immediate observations can be made that are quite 
general. First, ^sin may exceed ^com by several orders 
of magnitude. Second, in the case of weakly-reflective 
elements, it is necessary to use rather large values for 
N to achieve this orders-of-magnitude improvement in 
coupling strength. 



B. Linewidth narrowing 

The saturation phenomenon described above reveals 
another interesting effect as the number of elements is 
increased beyond A'opt : The presence of the array inside 
the cavity acts to narrow the cavity linewidth. The phys- 
ical basis behind this is rather transparent and relies on 
two observations. First, the fact that the array is trans- 
parent at the cavity resonance frequency means that the 
finesse of the cavity - which is related to the number of 
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FIG. 9. (a) Coupling strength, normalized to g, as a function 
of the number of elements. Four sets of data are shown. The 
red diamonds (blue circles) correspond to a single-element re- 
flectivity of 20% (99.4%). Closed (open) symbols represent an 
inter-element spacing of d = d_ (d = d- -I- 20 A). We note that 
the inter-element separation has a much stronger effect for 
highly-reflective elements. (L « 6.3 x 10*A, bare-cavity flnesse 
« 3 X 10*.) (b) A similar plot showing the effective cavity 
linewidth, normalized to the bare-cavity linewidth, obtained 
in each of the cases displayed in the upper plot, (c) Putting 
these two together we can calculate the cooperativity, normal- 
ized to the single-element cooperativity, and demonstrate an 
enhancement by several orders of magnitude for the chosen 
parameters. 



round-trips a photon makes inside the cavity on average 
- is unchanged by its presence. The second observation 
is that, as we have already noted, Eq. (38) reveals that 
the cavity is effectively lengthened to a length L^tf- Since 
the bare-cavity linewidth is Kc cc 1/{FL), it follows that 
the linewidth of the cavity is reduced to Keff oc 1/{J-Lcs). 
By choosing the right parameters, one can optimize for 
this linewidth-narrowing effect. Fig. 9(b), to narrow the 
optical resonance substantially. As several mechanisms 
in optomechanics, e.g., cooling in the linearized regime, 
are improved in the so-called 'resolved sideband' regime, 
where Wm ^ Kc, the collective mechanism we describe 
can result in the condition Kc > Wm S> n^s being satis- 
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FIG. 10. Reflectivity (red), transmission (green), and absorp- 
tion (orange) for N — 5 elements with an individual reflec- 
tivity of ca. 50% and an absorption of ca. 1.6% per element. 
The dashed blue curve is identical to the solid blue curve 
in Fig. 2. Note that the absorption around the transmissive 
points is highest close to d+ and lowest close to d_. 



fied, thereby improving the performance of these mecha- 
nisms in the system. 

One other figure-of-merit that is relevant to several 
mechanisms is the cooperativity C — g'^/{KT), i.e., 
the ratio of the square of the optomechanical coupling 
strength to the product of the optical and mechanical 
decay rates, k and F, respectively. For a multi-element 
system composed of independent oscillators and operat- 
ing in the transmissive regime, F is independent of iV, 
but (7sin ~ N^'^ and Kc ~ N^ in the appropriate regime; 
this is illustrated in Fig. 9(a) and (b). The result is a 
competition between these tvsro factors, yielding a con- 
stant cooperativity as N is increased. In Fig. 9(c) we 
plot the normalized cooperativity, i.e., the cooperativity 
for the A^-element ensemble divided by that for a single 
element; the enhancement obtained with the parameters 
used is of almost 10^. Even when absorption is included, 
as explained below in Sec. IV A, an enhancement of sev- 
eral orders of magnitude is still possible [58] . 



IV. RESILIENCE TO ABSORPTION AND 
IMPERFECTIONS 

In this section we address questions regarding the re- 
silience of the mechanism we introduce in the face of ex- 
perimental imperfections in the positioning, reflectivity, 
and absorption of individual elements. 



A. Eflfects of absorption in the elements 

For a nonzero per-element absorption, the largest 
amount of absorption appears close to the points where 
the ensemble is transparent. This is shown in Fig. 10, 
where we plot the reflectivity, transmission, and absorp- 
tion of an ensemble of 5 elements as the spacing between 



Im{C} = 
Im{C} = 10-= 
Im{C} = 10--' 
Iin{C} = 10-' 



6 8 10 12 14 

Number of membranes, A^ 



(b) 




6 8 10 12 14 

Number of membranes, A^ 

FIG. 11. Effect of increasing absorption on the linewidth 
(normalized to the bare-cavity linewidth) for a single-element 
reflectivity of 99.4% and (a) d = d-, (b) d = d-+ 20A. Larger 
inter-element separations make the system more tolerant to 
higher levels of absorption. The linewidth of the bare cavity 
is represented by the horizontal dashed black lines. These 
figures should be compared to the data shown in solid blue 
circles Fig. 9(b). (Bare-cavity finesse « 3 x lO"*, other param- 
eters as in Fig. 9.) 



the elements is scanned; this figure is meant to comple- 
ment Fig. 2. A general feature is that the absorption is 
largest at d_|_ and rather smaller at d_ ; we have therefore 
chosen the latter as our working point throughout this 
paper. 

A systematic study of the effect of absorption on the cav- 
ity linewidth is shown in Fig. 11. The cavity linewidth, 
both in the presence of absorption and in its absence, 
is calculated numerically by scanning over, and fitting a 
Lorentzian to, the cavity resonance. An approximate ex- 
pression for the linewidth in the presence of absorption 
can be given by taking into account the optical losses in 
the ensemble through a nonzero Im{C}, yielding 
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'i^off.abs 







where the factor A corresponds to the single-pass absorp- 
tion for the ensemble, and k^s is evaluated with C — >■ — |C|- 
For d = d-, Im{C} <^ 1, and TV > 2, ^ is given approx- 
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FIG. 12. Profiles of tlie 'sinusoidal mode' (i.e., the coupling 
strength of each element, in arbitrary units) in a 5-element 
array, for a representative sample of 1000 populations, for a 
total of 4000 populations, demonstrating random positioning 
errors of varying degree. The position of each element is var- 
ied from the ideal case by a random number drawn from a 
Gaussian distribution with standard deviation, from left to 
right, of IQ-^A, 10"^-^A, 10"^A, and lO'^'^A. The overall 
coupling strength resulting from each of these simulations is 
shown in Fig. 13. {( = —0.5, d = d-, other parameters as in 
Fig. 9.) 



imately by 



2 Ini{C} sin(i/) [Ccos(z^) + sin(i/)] Un_i [cos(7r/iV)] 



2iVIm{C} 



1 — cos( 



v/1 + lCpsinh 



N , 



COS 



arccos 



(f) 



arctan(ICI) 



arccotdCDKl)' (43) 



where v = kd- and U'-^{x) is the first derivative of 
the n*^ Chebyshev polynomial with respect to its argu- 
ment. Upon substituting this expression for A, Keff,abs 
agrees with the numerically-calculated data shown in 
Fig. 11. For large N and |C| (but Im{C} ^C 1) we find 
A^NlmiC}. 

To stay within the frame of the ID model considered here, 
a small misalignment in the individual elements can be 
modeled similarly to absorption, since both effects repre- 
sent a loss channel for the cavity field. Other detrimen- 
tal effects of absorption, such as heating, are mitigated 
by the large coupling strengths obtained, which allow 
much smaller photon numbers to be used {g^i^ oc N^ in- 
creases faster than the absorbed power as A^ increases). 
We note also that at large input powers it might be pos- 
sible to exploit photothermal forces to further enhance, 
or change the nature of, the collective optomechanical 
interaction [75-79]. 
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FIG. 13. Overall coupling strength for the populations shown 
in Fig. 12. Compared to the ideal case, the coupling strength 
for the four cases has an average (standard deviation) of 0.0% 
(1.0%), -0.3% (3.0%), -3.4% (9.4%), and -21.5% (37.2%). 
(Parameters as in Fig. 12.) 



B. Inhomogeneities in the system 

Here we discuss the behavior of the system in the face 
of two experimentally realistic sources of inhomogeneity: 
positioning errors and imperfections leading to different 
absorption in the elements making up the array. To ad- 
dress the issue of positioning errors, we performed nu- 
merical simulations with several thousand arrays {N — 5, 
C — —0.5, L « 6.3 X lO'^A, d = d-) whose elements were 
each shifted from the optimal position by a random shift 
drawn from a Gaussian distribution with a standard de- 
viation of 10 nm. The resulting coupling strengths were 
within ±10% of the analytically-calculated value. A rep- 
resentative sample of the analysis, extended to broader 
Gaussian distributions and performed on 4000 arrays, is 
shown in Fig. 12 and Fig. 13. 

Until now, for the most part, the elements making up the 
array were assumed to be identical and non-absorbing; 
the latter is an excellent approximation to a single mem- 
brane with an imaginary part of the refractive index be- 



ing 



< io~Mo- 



81]. Fluctuations in the polar- 



izability of the individual elements by up to ±10% give 
coupling strengths within ±6% of the analytic value. Our 
numerical investigations show that the major effect of ab- 
sorption is not on the optomechanical coupling strength, 
but on the linewidth of the cavity. High levels of absorp- 
tion would therefore not alter the coupling strength sig- 
nificantly, but would limit the achievable cavity finesse. 



V. NUMERICAL EXAMPLES 

In this section we give two numerical examples, one us- 
ing readily-available membranes and the other state-of- 
the-art highly-reflective ones, to demonstrate the power 
of transmissive optomechanics. If we take commercial sil- 
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icon nitride membranes [72, 80, 82, 83] with an intensity 
reflectivity of 20% (C = —0.5) and Xq = 1.8 fm, and a 
cavity with L — 6.7 cm and a wavelength of 1064 nm, we 
can estimate gcom « 27r x (12.8 x N'^^^ Hz) for N > 3. 
For the sinusoidal mode, and with the same parameters, 
gsin ~ 27r X (1.3 X N^/^ Hz) for large N; an improvement 
by over an order of magnitude when A^ = 10 (cf. Fig. 7). 
Suppose, now, that we consider instead highly-reflective 
membranes [73, 74] having 99.4% intensity reflectivity 
(C = -12.9), xo = 2.7 im, and Wm = 27r x 211kHz. For a 
0.25 cm-long cavity with finesse F = 1.2 x 10^, d = d^, 
and N — 5 membranes, one obtains (/com ~ 27r x 600 Hz 
and gsin ~ 27r x 270 kHz, which is larger than both w^ 
and Kc = 27r X 250kHz. At a temperature of IK and 
with a mechanical quality factor of 10^ the single-photon 
cooperativity is ca. 14 for this system; strong coupling 
between a single photon and a single phonon is already 
within reach with only a few elements. 

Moving away from highly- reflective scatterers, we can 
apply our results to systems of very low reflectivity, such 
as atoms, molecules, dielectric microspheres, etc. It 
should be first noted that all these systems have a reflec- 
tivity on the order of 10~^, which means that N\(^\ ^ 1 in 
typical experiments (for example with cold atoms in cav- 
ities [67-70]), i.e., the particles do not significantly mod- 
ify the mode structure of the cavity resonance. In this 
case gsin reduces to the expected N^^"^ scaling that arises 
from the independent coupling of well-localized scatter- 
ers interacting with an unmodified cavity field [67-70]. 
We note, however, that recent experiments [84] using 
"pancake" -shaped clouds of cold atoms in an optical lat- 
tice have shown intensity refiectivities as high as 80% 
and are approaching a regime where the effects discussed 
previously may be observed. 



VI. CONCLUDING REMARKS 

Transmissive optomechanics presents a departure from 
traditional optomechanical systems in that the refiectiv- 
ity of a compound element is purposely engineered to be 
as close to zero as possible. A wealth of interesting ef- 
fects exist in this regime, not least (i) the possibility of 
strongly increasing the optomechanical cooperativity and 
obtaining strong coupling between a single photon and a 
single phonon, (ii) the existence of a linewidth-narrowing 
mechanism that renders the resolved-sideband regime of 
optomechanics more accessible, (iii) the existence of long- 
range interactions within optomechanical arrays [58] , and 
(iv) the possibility to enhance optomechanical nonlinear- 
ities [40]. The system we described may be composed of 
any periodic array of linearly- interacting polarizable scat- 
terers, e.g, an ensemble of macroscopic dielectric scatter- 
ers, or even atoms in an optical lattice, and therefore 
presents a widely-configurable and robust basis for inves- 
tigating strong and collective effects in optomechanical 
or electromechanical setups. 
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FIG. 14. Modeling an individual membrane as a plate of 
thickness I and refractive index n. 
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Appendix A: Polarizability of thin dielectric plates 

We model each membrane, or any other thin dielec- 
tric scatterer, as a plate of (real) refractive index n and 
thickness I (see Fig. 14). Going from left to right, the 
two interfaces have amplitude reflectivities ti = r = 
(1 — n)/{l + n) and r2 = — r, as well as transmissivi- 
ties ti — 2/(1 + n) and t2 = 2n/(l -I- n), which are ob- 
tained from the respective Fresnel coefficients at normal 
incidence. The amplitude reflection and transmission co- 
efficients for light of wavenumber k = uj/c incident on 
the plate as a whole can then be written 



r 1 



^ — 2inkl\ 



1 — r^e 



-2inkl 
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(Al) 



respectively. We then define the effective membrane po- 
larizability ^ = —ir/t [59], whereupon 



2r 

C — siiiinkl) 



2n 



■ sin^nkl) . 



(A2) 



The 'thin membrane' approximation corresponds to tak- 
ing n ^ 1 but nkl <?C 1, whereby we may approximate 



c 



-^n^kl. 



(A3) 



This justifies our assertion throughout this paper that 
C < for lossless scatterers. Let us now calculate the 
relative change in C over a frequency interval of size duj — 
cSk centred around the frequency uj: 
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Throughout the main text, we were interested in the op- 
tomechanical behavior of arrays of such structures inside 
cavities. We may then reasonably take 5uj/uj = l/Q, 
where the quality factor Q of an optical cavity is gener- 
ally of the order of 10^-10^. Under the thin membrane 
approximation one can write 



dC/dk 



5k 



1 
Q 



(A5) 



For membranes that are thin on the scale of a wavelength, 
this requirement is thus satisfied trivially; for thicker 
membranes one must take care to avoid the structural 
resonances at sin(nfcl) = 0. Another reason for requiring 
the membranes to be thin is that the free-spectral range 
of each 'membrane-etalon' must be significantly larger 
than that of the main cavity, which requires I <SS^ L/n, 
where L is the length of the cavity. A final, and stronger, 
requirement is that the entire membrane stack should fit 
within the Rayleigh range of the cavity mode. 



where 7 is the (half-width at half-maximum) linewidth 
of the atomic resonance and A = wl — Wa the detun- 
ing between the driving field (frequency wl) and atomic 
resonance (i^Ja)- o'a, — 3Aa/(27r) is the resonant radiative 
cross-section of the atom at a wavelength Aa — 27rc/a;a. 
A first-principles derivation of this relation can be found 
in Ref. [64]. A 'non-absorbing' scatterer, where C. is real, 
can be approximated by a TLA when |A| ^7. 
The expression above assumes that one is working in the 
regime of low saturation, where the excited state popu- 
lation is negligible. Under these conditions, the effect of 
having m atoms inside the waist of the Gaussian beam 
is simply ^ — >■ m(. 



Appendix B: Polarizability of two-level atoms 



For a single two-level atom (TLA) at the waist of a 
Gaussian field of mode area ctl, we may write 



c = - 



7 



2o-L iA + 7 ' 



(Bl) 
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